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We construct a spin- drift- diffusion model to describe spin-polarized electron transport in 
zincblende semiconductors in the presence of magnetic fields, electric fields, and off-diagonal strain. 
We present predictions of the model for geometries that correspond to optical spin injection from the 
absorption of circularly polarized light, and for geometries that correspond to electrical spin injection 
from ferromagnetic contacts. Starting with the Keldysh Green's function description for a system 
driven out of equilibrium, we construct a semiclassical kinetic theory of electron spin transport in 
strained semiconductors in the presence of electric and magnetic fields. From this kinetic theory 
we derive spin-drift-diffusion equations for the components of the spin density matrix for the spe- 
cific case of spatially uniform fields and uniform electron density. We solve the spin-drift-diffusion 
equations numerically and compare the resulting images with scanning Kerr microscopy data of 
spin-polarized conduction electrons flowing laterally in bulk epilayers of n-type GaAs. The spin- 
drift-diffusion model accurately describes the experimental observations. We contrast the properties 
of electron spin precession resulting from magnetic and strain fields. Spin-strain coupling depends 
linearly on electron wave vector and spin-magnetic field coupUng is independent of electron wave 
vector. As a result, spatial coherence of precessing spin flows is better maintained with strain than 
with magnetic fields, and the spatial period of spin precession is independent of the applied electrical 
bias in strained structures whereas it is strongly bias dependent for the case of applied magnetic 
fields. 

PACS numbers: 72.25.Dc, 71.70.Ej, 85.75.-d 



I. INTRODUCTION 



A new generation of electronic devices, with the po- 
tential to outperform conventional electronic circuits in 
speed, integration density and power consumption, has 
been proposed based on the ability to manipulate elec- 
tron spin in semiconductors [1-3]. To design semicon- 
ductor structures whose function is based on electron 
spin, it is necessary to understand spin dynamics and 
spin-polarized transport, and in particular, how they are 
affected by electric, magnetic and strain fields. Spin 
dynamics and spin transport in semiconductors have 
been studied experimentally using time- and/or spatially- 
resolved spin-sensitive optical spectroscopies based on 
the magneto-optical Faraday and Kerr effects [4-17]. 
Theoretical approaches to describe spin dynamics and 
spin-polarized transport include: the two-component 
drift-diffusion model [18-21], Boltzmann equation ap- 
proaches [22-31], Monte-Carlo techniques [32-34] and 
microscopic approaches [35-37]. Many of these works 
have focused on the description of spin-relaxation times. 
In this paper we are particularly interested in describing 
the effects of strain on spin-polarized electron transport. 

Recently, the possibility of using strain to control elec- 
tron spin precession in zincblende structure semiconduc- 
tors has been demonstrated [13-17]. The spatial part 
of the conduction electron wave function is modified in 
the presence of strain which affects the electron spin de- 
grees of freedom due to spin-orbit coupling. Scanning 



polar Kerr-rotation microscopy (in which a linearly polar- 
ized laser propagating along the sample normal is raster 
scanned in the x-y sample plane and the polarization ro- 
tation angle of the reflected beam is detected) has been 
shown to be a valuable experimental technique to image 
electron spin flow in zincblende structure semiconduc- 
tors. Scanning Kerr microscopy has recently been used to 
acquire 2-D images of electron spin flows in n-type GaAs 
epilayers subject to uniform electric, magnetic, and (off- 
diagonal) strain fields in which spin-polarized electrons 
have been optically injected by circularly polarized light 
[15] or electrically injected into lateral spin-transport de- 
vices using ferromagnetic iron (Fe) contacts [16]. 

In this paper, we microscopically derive an equation 
of motion for the electron Green's function giving a full 
quantum-mechanical description of electron spin dynam- 
ics and transport in the presence of electric, magnetic and 
strain fields. From this equation of motion we construct 
a semiclassical kinetic theory of electron spin dynamics 
and transport in the presence of these fields. From the 
semiclassical kinetic theory, we derive a set of spin-drift- 
diffusion equations for the components of the spin den- 
sity matrix for the case of spatially uniform fields and 
electron density. This approach extends beyond the two- 
component drift-diffusion model, because the spin distri- 
bution function has the form of a spin-density matrix to 
account for spin coherence effects. We are particularly in- 
terested in describing spin polarized electron transport in 
the presence of momentum-dependent coupling between 
spin and off-diagonal strain. Off-diagonal strain arises, 
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e.g., from uniaxial stress along a (110) crystal axis. Be- 
cause the spin-strain coupling depends on electron wave 
vector, the diffusion coefficient and mobility appear in the 
strain coupling terms of the spin-drift-diffusion equation. 
We solve the spin-drift-diffusion equations numerically 
and find good agreement with scanning Kerr microscopy 
images of spin-polarized conduction electrons flowing lat- 
erally in bulk epilayers of n-type GaAs. We contrast 
the effects of magnetic and strain fields on electron spin 
transport. Both magnetic fields and strain fields (with 
off-diagonal strain components) lead to electron spin pre- 
cession in zincblende semiconductor structures. How- 
ever, spin- strain coupling is linear in electron wave vec- 
tor whereas spin-magnetic field coupling is independent 
of electron wave vector, leading to qualitatively differ- 
ent behavior. We present calculations that contrast spin 
transport in strain and magnetic fields and then present 
a series of theoretical predictions for (and calculated im- 
ages of) spin transport in geometries that correspond to 
optical spin injection and in geometries that correspond 
to electrical spin injection into lateral devices. 

The outline of the paper is as follows: in Sec. II, we 
solve the spin-drift-diffusion equations for the spin den- 
sity matrix (derived in the Appendix) and we compare 
results with experimental data; in Sec. Ill, we discuss 
differences between strain-induced and magnetic-field- 
induced spin precession; a series of theoretical predic- 
tions for spin transport are presented in Sec. IV; and 
in Sec. V, we summarize our conclusions. Details of the 
derivation of the transport equations are presented in the 
Appendix. 

II. SPIN-DRIFT-DIFFUSION EQUATIONS 

Figure 1 shows schematics of the two experimental ge- 
ometries considered. A [001] -oriented n:GaAs epilayer 
sample, whose surface normal is along the z-axis, is sub- 
ject to in-plane electric and magnetic fields, as well as 
strain fields. As illustrated in Fig. la, the sample may 
be optically excited by a circularly polarized laser beam 
propagating along the surface normal so that electrons, 
spin polarized along the z-axis, are optically injected. Al- 
ternatively, spin polarized electrons may be electrically 
injected into the semiconductor from ferromagnetic (FM) 
contacts as illustrated in Fig. lb. For electrical injec- 
tion, the injected electron spin polarization follows the 
magnetization M of the ferromagnetic contact, and is 
therefore typically in the x-y plane of the sample surface 
(for certain special ferromagnetic contacts both M and 
the injected electron spin polarization may be directed 
along the surface normal). The spin polarized electrons 
subsequently drift and diffuse in the x-y sample plane. 
The epilayer, in which the electrons are confined, is thin 
compared to a spin diffusion length so that the spin den- 
sity is essentially uniform in the z-direction. The result- 



ing steady state spin polarization is imaged via scanning 
Kerr microscopy, which is sensitive to the z-component of 
electron spin polarization, n^. In the experiments [15,16], 
ohmic side contacts allow for a lateral electrical bias {Ex) 
in the x-y plane, Helmholtz coils provide an in-plane mag- 
netic field {By)^ and uniform off-diagonal strain (cxy) is 
provided by a controlled uniaxial stress applied to the 
sample along a [110] crystal axis. The experimental data 
in the lower panels of Fig. 1 show false color 2-D im- 
ages of the measured z-component of electron spin polar- 
ization n^, showing the flow of spin polarized electrons. 
In the lower left panel, spin-polarized electrons are opti- 
cally injected using a laser focused to a 4 /im spot (spot 
size illustrated by the circular red dot) and these spins 
subsequently diffuse and drift to the right in an applied 
electric field {E^ = 10 V/cm). In the lower right panel 
spin-polarized electrons are electrically injected from an 
iron (Fe) tunnel-barrier contact that is magnetized along 
—X. The injected spin polarization is in the x-direction, 
and a small magnetic field in the y-direction {By = 3.6 
G) rotates the spins so that there is a z-component of spin 
polarization that can be detected by Kerr microscopy. 
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FIG. 1. Diagrams and corresponding experimental data of 
the two measurement geometries considered in this work, (a) 
Spin polarized electrons, oriented along are optically in- 
jected into an n:GaAs epilayer with a focused, circularly po- 
larized laser incident along the z-axis (indicated by the red 
arrow). Lateral contacts provide an in-plane electric field 
(E II ±x), Helmholtz coils provide an in-plane magnetic field 
(B II and off-diagonal strain Cxy results from controlled 

uniaxial stress applied along the [110] GaAs crystal axis. 
Scanning polar Kerr-rotation microscopy is used to spatially 
image Uz , the z-component of the resulting steady-state elec- 
tron spin polarization. The drift, diffusion and precession of 
the spin flow are determined (in part) by the applied elec- 
tric, magnetic and strain fields (see data below; adapted from 
Ref. [15]). (b) Electrical injection of spin-polarized electrons 
into an niGaAs epilayer using ferromagnetic (FM) injector 
contacts. The contact magnetization (M) may be either in 
the x-y plane, or out-of-plane (M||z). As shown schemati- 
cally, the injected spins then flow laterally within the niGaAs 
epilayer, and riz can be spatially imaged with scanning Kerr 
microscopy. Experimental data, using an iron injector with 



M||- 
[16]). 



X and By = 3.6 G is shown below (adapted from Ref. 



To theoretically describe the effects of spin-strain cou- 
pling, we consider a Hamiltonian of the form 

H=^[p- f A(R)]2 + em) + f ^^^^ • ^, (1) 

(Details of our microscopic treatment, that also include 
Dresselhaus [38] and Rashba [39] spin coupling terms, 
are described in the Appendix.) Here a is the vector 
having the three Pauli spin matrices as its components, 
e = — |e| is the charge of an electron, m is the conduction 
electron eflFective mass, (j) is the scalar potential, g is the 
conduction electron gyromagnetic ratio, /is is the Bohr 
magneton, B is the applied magnetic field and A is the 
corresponding vector potential. The term 



(2) 



where (fy and ifz are obtained by cyclic permutations, 
describes the lowest-order k-linear coupling of electron 
spin to the off-diagonal components of the crystal strain 
tensor [40]. Here, eij is the strain tensor, 



-h 



C2V 



S'^^2Egm{l-rj/3) ' Eg + A 



C2 is the interband deformation potential constant for 
acoustic phonons, A is the spin-orbit splitting of the va- 
lence band and Eg is the band gap energy. For electrons 
moving laterally in the x-y sample plane, the spin-strain 
coupling term takes the form exy{(Jxky — cFykx)^ which has 
the same form as the Rashba spin coupling in asymmet- 
ric heterostructures. That is, the presence of oflF-diagonal 
strain results in an eflPective magnetic field Bg which is 
in-plane, orthogonal to the electron momentum k, and 



grows linearly with the magnitude of k. In contrast, spin- 
magnetic field coupling is independent of k. Because of 
this diflFerence, strain and magnetic fields aflFect the flow 
of spin-polarized electrons in qualitatively diflFerent ways. 

We derive spin- drift- diffusion equations for the spin 
density matrix 



Ps = n(R) + ^ niilVjcTi 



(3) 



=x,y,z 



accounting for the quantum nature of spin. Here n is the 
total electron density and rii denote the three spin density 
components. Information about spin polarization along 
any given direction s is found by Ti ps "^i^x y z ^i^i- We 
consider a case, corresponding to the experimental ge- 
ometry, in which the net electron density is position in- 
dependent and the electric, magnetic and strain fields 
are uniform. As shown in the Appendix, for electric and 
magnetic fields in the x-y plane and stress along a [110] 
crystal axis, the system of spin-drift-diflFusion equations 
for the spin density matrix is given by 



DV^ -h /iE • Vr - C'^D 



1 



DV'^ -h /iE • Vr - ClD 



riz 



-G.(R), (4) 



Cs ( 2D^+pEy 



-G,(R), (5) 



- 2ClD - -\ 


















Uy = 



ny = -G,(R), (6) 



where is the spin relaxation time, D is the electron 
diflFusion coefficient, p is the electron mobility, Cs = 
Csmcxy/h'^^ hCsi — Ql^sBi {i = x^y) and the vector 
G describes generation of the spin components. Because 
the spin-strain coupling depends on electron wave vec- 
tor, both D and p appear in the strain coupling terms. 
D and p do not appear in the magnetic field coupling 
terms. A second order term in strain, C^D, appears in 
the drift diffusion equations. The analogous second order 
term in magnetic field is proportional to the momentum 
relaxation time, which is short on the time scales of in- 
terest here, and thus this term is very small and can be 
neglected (see the Appendix for details.) 

Performing a 2-D Fourier transform to momentum 
space, the system of equations for the spin polarization 
components, in our experimental geometry, is given by 
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Oi(q)nx(q) + 02(q)n,(q) = -G,(q) , (7) 
Oi(q)n„(q) + 03(q)n.(q) = -Gyiq) , (8) 
04(q)n.(q) - 02(q)n,(q) - 03(q)%(q) = -G,(q), (9) 



where 



Ol(q) 



-Dq^ + i/iE • q - C^D 



02(q) = [Cb, + a(i2Z)fe + M^x)] , 
03(q) = -[Cb^ - Cs{i2Dqy + nEy)] , 



04 (q) 



-Dq^ + i/iE ■ q - 2C,D 



(10) 

(11) 
(12) 

(13) 



The resulting spin density, n^(R), is then given by the 
2-D Fourier transform of ^^(q). 



i(R) = ^ y e'^-^n,iq)d\ 



(14) 



The solutions for the spin polarization components 
from this system of equations are 

, , -(Q4 + ^)Gx + (^)G, + 02G. 
nx(q) = , (15) 



(O4 + %^)Gy + O3G, 



nziq) 



02Gx + 03Gy + OiG, 
d 



where the denominator, d, is 

d = 0,04 + (O2)' + (Os)^. 



(16) 



(17) 



(18) 



In Fig. 2 we compare our calculations with experimen- 
tal results for the case of optical spin injection. These 
experimental data were obtained via scanning Kerr mi- 
croscopy [15]. In these measurements, a linearly polar- 
ized narrowband Ti:sapphire laser, tuned just below the 
GaAs band-edge and focused to a 4 /im spot on the sam- 
ple surface, was raster-scanned in the x-y sample plane 
to construct a 2-D image of n^, the z-component of the 
electron spin density. The sample was a Si-doped n-type 
GaAs epilayer (electron density = 1 x 10^^ cm~^) grown 
on a [001] oriented semi-insulating GaAs substrate. Spin 
polarized electrons were optically injected into the sam- 
ple by a separate, circularly-polarized 1.58 eV diode laser 
that was also focused to a 4 /im spot on the sample. Mea- 
surements were performed at a temperature of 4 K. For 
additional experimental details, see Ref. [15]. Figures (a), 
(c), and (e) compare calculated and measured results for 
spin flowing to the right in the presence of uniform off- 
diagonal strain {cxy = 4 x 10~^; Ex = 12 V/cm), while 
Figs, (b), (d), and (f) compare calculated and measured 
results for spins flowing to the right in the presence of an 



applied magnetic field {By = 16 G] = 7 V/cm). There 
is very good agreement between theory and experiment. 
The material parameters used in the calculation were: 
mobility /i = 3000 cm^/Vs, spin relaxation time = 125 
ns (we use these values for /i and throughout the pa- 
per), and diffusion constant D = 10 cm^/s. Throughout 
the paper we use a value of Cs = 4.0 eVA for the spin- 
strain coupling coefficient based on the experiments of 
Ref. [41]. Both strain and magnetic field lead to preces- 
sion of the electron spins. However, the spatial damp- 
ing of the precession is more pronounced when magnetic 
rather than strain fields are applied (at the same preces- 
sion length). In the remainder of the paper, we discuss 
various predictions of the spin-drift-diffusion model. 
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FIG. 2. (a) and (b) show images of the calculated 
z-component of electron spin density (riz) for optically in- 
jected spins that are flowing to the right and that are subject 
to (a) applied [110] stress (strain Cxy = 4 x 10~^; Ex = 12 
V/cm), and (b) applied magnetic field {By = 16 G; Ex = 7 
V/cm). Note that the color scale is modified such that black 
indicates negative values of riz (i.e., spins that have precessed 
into the sample plane). For all calculations, D = 10 cm^/s, 

= 3000 cm^/Vs, and Ts = 125 ns. (c) and (d) show images 
of corresponding experimental data under the same condi- 
tions, (e) and (f) directly compare line cuts through the cen- 
ter of the corresponding calculated (black line) and measured 
(red dots) images. 
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III. CONTRASTING THE EFFECTS OF 
MAGNETIC AND STRAIN FIELDS 



Although strain and magnetic fields both lead to elec- 
tron spin precession, spin precession due to strain dif- 
fers qualitatively from that due to a magnetic field be- 
cause spin-strain coupling depends linearly on the elec- 
tron wave vector k whereas spin-magnetic field coupling 
is independent of electron wave vector. For spin preces- 
sion in an applied magnetic field there is a pronounced 
spatial dephasing of spin precession due to the random- 
izing nature of diffusion. The net spin polarization at 
a remote location from the point of spin generation re- 
sults from the combined sum of many random walk paths. 
Each path takes a different amount of time, giving a dif- 
ferent precession angle. As a result the spatial coherence 
of the spin flow is rapidly lost. By contrast, the spin 
precession frequency in a strain field scales linearly with 
the electron wave vector. As a result the precession an- 
gle is correlated with electron velocity and therefore with 
the electron's position. There is a partial cancelation of 
accumulated precession angle during the parts of the dif- 
fusive path that are traversed in the opposite direction, 
due to the linear dependence of effective magnetic field 
on electron momentum. Indeed, for motion strictly in 
one dimension, the spin ensemble would not dephase at 
all (however, it would still decohere with a timescale Tg). 
Scattering from +k to — k simply reverses the direction of 
precession, leading to an exact correspondence between 
spatial position and spin orientation. 

A second difference between spin precession due to 
strain and that due to magnetic field is that for spin- 
magnetic field coupling the spatial precession length 
increases with increasing electric field (drift velocity), 
whereas for spin-strain coupling this length scale is inde- 
pendent of electric field. For spin-magnetic field coupling, 
the spin precession frequency is independent of electric 
field, but because the electrons drift faster in a larger 
electric field, the spatial precession length increases. For 
spin- strain coupling, the precession frequency scales di- 
rectly with the average electron velocity, so that the spa- 
tial spin precession length is nearly independent of elec- 
tric field. 

To illustrate these differences, we consider the sim- 
ple geometrical case in which spins are generated by a 
stripe optical excitation along the y-axis, the electric field 
E II X, the magnetic field B || ^ and off-diagonal strain 
exy is generated by stress applied along x (a [110] crys- 
tal axis). For this geometry, spin-drift-diffusion occurs 
along the x-axis so that the spin densities do not de- 
pend on the y-coordinate. Because of the orientation of 
strain and magnetic field, Uy vanishes. For this geomet- 
rical case, O3 vanishes and only the x-component, is 
of interest. The Green's function for the z-component 
of spin density precesses and decays exponentially as 



9z{x) = Aexp(7^ ^^^ ) -\- CC. The four complex decay 
parameters, ry, can be found from the four zeros of the 
denominator d in Eq. (18), 77 = iq^y/Dr, 
a zero of d. Setting a = (iExTs/ \/Dt. 
ps = CsVDFs, Q 
of d are given by 



s, where Qx is 



(|)^ + 1, and C = 7y + f , the zeros 



a 
2 • 



s •) 



[e -{Q + [C' - (Q + Ws^)] = -[-/3b + WsCf 



(19) 



For strictly 1-D motion, in which the electrons are con- 
strained to move only along the x-axis, the spin-drift- 
diffusion equations are modified in that the factor of 2 
multiplying Cs^D in Eq. (6) is replaced by unity [see 
Eq. (48)]. The corresponding factor of 2 that appears in 
04(q), Eq. (13), and that multiplies (3s^ in the second 
factor on the left hand side of Eq. (19) are also replaced 
by unity for this strictly 1-D case. The magnetic field 
terms of the spin- drift- diffusion equations are the same 
for the 1-D and 2-D cases. 

We consider the solutions of Eq. (19) for three cases: 
(1) applied magnetic field without strain; (2) strain with- 
out magnetic field; and (3) strain without magnetic field 
for strictly 1-D motion. For case (1) the solutions are. 



(20) 



For case (2) the solutions are. 



c = ±i 



VQ + /?.V0 + 2/3." + (Q-^) 



VQ + Pe' VQ + W-iQ-^) ^21) 



For case (3) the solutions are, 



(22) 



From Eq. (20) we see that the real part of ^ increases 
with increasing magnetic field so that the decay length 
for spin density decreases with increasing magnetic field. 
From Eq. (21), we see that the real part of ( also in- 
creases with increasing strain, but that there is a partial 
cancellation, due to the appearance of the term — 
The decay length for spin polarization also decreases with 
increasing strain, but more slowly than for magnetic field. 
To compare the strain and magnetic field cases, we con- 
sider strains and magnetic fields that give the same spa- 
tial precession length. For the case of strain in strictly 
1-D motion, the factor of 2 and of 1/2 that appear in the 
numerator of Eq. (21) are both replaced by unity and the 
partial cancelations that occur for the 2-D case become 
complete in the 1-D case, giving the result of Eq. (22). 
In this case strain does not reduce the decay length of 
spin polarization. 
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In Fig. 3, we calculate as a function of position for 
a narrow stripe excitation. The solid curve, which is the 
same in each of the panels, is for the non-precessing case 
of electric field only, with neither strain nor magnetic 
field. The dashed lines in the upper row show the case 
of applied magnetic field. The dashed lines in the cen- 
ter row show the case of applied strain, and the dashed 
lines in the lower row show the case of applied strain for 
strictly 1-D motion. If there is no strain, Eqs. (4)- (6) 
for the strictly 1-D motion case are the same as the 2-D 
case. The values of magnetic field and strain were chosen 
so that the spatial precession periods are approximately 
equal for the two cases. The values of strain and mag- 
netic field in the right column of Fig. 3 are twice that in 
the left column. Comparing the solid and dashed lines 
in the upper row of Fig. 3 shows that decreases sig- 
nificantly when a magnetic field is applied. When the 
magnetic field is increased, the decay length is shorter. 
When strain is applied instead of magnetic field, as in 
the center row of Fig. 3, Uz also decays more rapidly. 
However, the decay length remains longer for the case of 
strain than for the case of an equivalent magnetic field. 
When strain is applied for the case of strictly 1-D mo- 
tion, the decay length is unchanged. The accumulated 
spin-precession angle due to strain depends only on the 
electron's location and not on the path taken to reach 
that location, resulting in no dephasing of the ensem- 
ble spin-polarization. This behavior does not occur for 
strictly 1-D motion with magnetic field. 




X (i^m) 

FIG. 3. Calculated z-component of spin polarization (uz) 
generated by a narrow stripe excitation at x = under the 
influence of a flxed electric field of Ex = 15 V/cm and: Top 
row, a magnetic field of = 15 G (left) or By = 30 G 
(right); center and lower rows, strain Cxy = 3.5 x 10~^ (left) 
or Cxy — 7 X 10~^ (right). The solid line is for electric field 
alone, without magnetic field or strain. The dashed lines rep- 
resent cases where, in addition to the electric field, a magnetic 
or strain field is applied. The bottom row is for the case of 
strain and strictly 1-D motion. D = 10 cm^/s, ji = 3000 
cm^/Vs, and Ts = 125 ns. 

In Fig. 4, we compare the eff"ect of increasing the 
electric field while in the presence of either a fixed mag- 
netic field or a fixed off-diagonal strain. The geometry 
and material parameters are the same as for Fig. 3. 
The two panels in the left column show calculated re- 
sults for 5^ = 15 G and Ex = 10 V/cm (top panel) 
and Ex = 20 V/cm (bottom panel). The two panels in 
the right column show calculated results for a strain of 
^xy = 3.5 X 10~^ and the same electric fields as in the 
left column. Comparing the top and bottom panels in 
the left column, we see that for fixed applied magnetic 
field the spatial precession period increases with increas- 
ing electric field. By contrast, comparing the top and 
bottom panels in the right column, we see that for fixed 
strain the spatial precession period does not change with 
increasing electric field. This behavior was observed in 
Ref. [15]. For both strain and magnetic field, the spatial 
decay length increases as the electric field increases, due 
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to the fact that spin transport is governed more by drift 
than by diffusion. 
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compare and contrast spin flows in the case of smah dif- 
fusion constant (left column; D = 1 cm^/s) with the case 
of large diffusion constant (right column; D = 20 cm^/s). 
In Fig. 5(a) and (b), there is no applied magnetic fleld 
and no effective magnetic fleld due to strain, therefore 
the spin density decays monotonically away from the 
point of injection. The larger diffusion constant in Fig. 
5(b) is reflected in the greater diffusive spread of spin po- 
larization along the y-axis. In Figs. 5(c) and (d) there is 
an additional magnetic fleld {By = 20 G), and the spins 
now precess as they flow to the right. With increasing 
distance from the point of spin injection, many spin pre- 
cession cycles are observed when the diffusion constant is 
small (Fig. 5c). In contrast, only a single full precession 
cycle is visible when the diffusion constant is large (Fig. 
5d). These images directly reveal how the effects of an 
increased diffusion constant can lead to signiflcant spatial 
dephasing of precessing spin flows. Figures 5(e) and (f) 
show the spin flow precessing in the presence of an effec- 
tive magnetic fleld due to strain {cxy = 7 x 10~^). Again, 
the degree of spatial spin coherence is much higher when 
the diffusion constant is smaller. Further, the line-cuts in 
Figs. 5 (g) and (h) show that the degree of spatial spin 
coherence is better maintained for the case of strain than 
for the case of magnetic fleld, as discussed in section III. 



D=1 cm2/s 



D=20 cm2/s 



FIG. 4. Calculated z-component of spin polarization (uz) 
generated by a narrow stripe excitation at x = under the 
influence of a lateral electric field Ex = 10 V/cm (top row) 
and Ex = 20 V/cm (bottom row). = 15 G in the left 



column: 



3.5 X 10 in the right column. The solid line 



is for electric field alone, without magnetic field or strain (no 
precession). The dashed lines represent cases where, in addi- 
tion to the electric field, a magnetic field or strain is present. 
D = 10 cmVs, M 3000 cm^/Vs, and Ts 125 ns. 



IV. PREDICTIONS OF THE 
SPIN-DRIFT-DIFFUSION MODEL 

In this section we present specific examples and predic- 
tions of the spin-drift-diffusion model. We flrst consider 
geometries that correspond to optical spin injection (ini- 
tial spin density no || ^), and then consider geometries 
that correspond to electrical spin injection (no || with 
spatially extended source contacts). 

Figure 5 shows the deleterious effect of diffusion on the 
spatial coherence of precessing spin flows. A steady-state 
spin polarization no || i is generated by a Gaussian source 
term with 4 /im full- width at half-maximum (FWHM), 
similar to the case of local optical spin injection. Spin- 
polarized electrons diffuse, and drift to the right in the 
presence of a lateral electric fleld, = 10 V/cm. We 
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FIG. 5. Calculated 80 x 120 /im images of Uz for opti- 
cally-injected spins, showing how diffusion directly influences 
the degree of spatial spin coherence for precessing spin flows. 
Images in the left column are for D = 1 cm^/s; images in 
the right column are for D = 20 cm^/s. All calculations use 
^ = 10 V/cm, /i 3000 cm^/Vs, and Ts 125 ns. (a) 
and (b) show spin flow in the absence of magnetic or strain 
fields, (c) and (d) show precessing spin flow in the presence 
of By = 20 G. (e) and (f) show precessing spin flow in the 
presence of strain (cxy = 7x 10~^). (g) and (h) show line-cuts 
through the center of the corresponding images, with blue and 
red lines showing the case of finite By and e^y, respectively. 



Figure 6 shows how the effective magnetic field due 
to strain (Be) can either augment or oppose a real mag- 
netic field (By), leading to markedly different spin flows 
depending on the flow direction. The schematic diagrams 
at the top of Fig. 6 illustrate how a uniform applied mag- 
netic field By is vectorially summed with the effective 
magnetic field Be due to finite strain, for the case of spin 
polarized electrons diffusing radially away from the point 
of injection (no drift). Since Be is always orthogonal to 
electron momentum k, Be is chiral for radially-diffusing 
spins, as shown. Under the conditions of Fig. 6, the 
vector sum of these two fields is given by B^et, which 
is finite for spins diffusing to the left, and negligible for 
spins diffusing to the right. In Fig. 6(a), spin polariza- 
tion no II i is generated by a 4 /im Gaussian FWHM 
source term, the applied field By = 6 and the strain 
is = 2 X 10~^. Spin polarized electrons are subject 
to diffusion only {Ex = 0). The resulting image of Uz is 
completely asymmetric - spins diffusing to the left "see" 
a net magnetic field and precess, while electrons diffusing 
to the right experience negligible net magnetic field and 
do not precess. This behavior was observed in Ref. [15]. 
The asymmetry is even more apparent when the spins 
drift in the presence of a lateral electric field E || ±x. 
Fig. 6 (b) shows that spins flowing to the left {Ex = — 10 
V/cm) exhibit pronounced precession with many preces- 
sion cycles, while Fig. 6 (c) shows that spins flowing 
to the right under the opposite lateral bias {Ex = +10 
V/cm) do not precess {By and Be effectively cancel in 
this direction, under these conditions). Fig. 6 (d) shows 
corresponding line-cuts through the images. 
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FIG. 6. Diagram and calculated images showing how an 
effective magnetic field due to off-diagonal strain (Be) can 
either augment or oppose a real magnetic field By, depend- 
ing on the electron momentum k. (a) Calculated 60 x 60 
/im image of spin diffusion with no electric field, and with 
By = 6 G and strain Cxy — 2 x 10~^. Spins diffuse radially 
away from the point of injection, "sampling" all momenta k 
in the x-y plane. Spins diffusing to the left see a net effec- 
tive magnetic field and precess; spins diffusing to the right 
see negligible net effective magnetic field and do not precess. 
(b) With Ex = —10 V/cm, electrons flowing to the left show 
pronounced spatial precession, (c) With Ex = +10 V/cm, 
electrons flowing to the right show essentially no precession, 
(d) Corresponding line-cuts through the center of the images 
in (b) and (c). In all images, D = 3 cm^/s, Ts = 125 ns, and 
fi 3000 cmVVs. 

The spin drift-diffusion equations [Eqs. (15) - (18)] 
can also be used to compute the in-plane components of 
electron spin polarization {ux and Uy). Knowledge of Ux 
and Uy is useful for comparison with experiments that 
are directly sensitive to in-plane components of electron 
spins, such as longitudinal or transverse (as opposed to 
polar) magneto-optical Kerr effect studies, or for direct 
electrical measurement of in-plane electron spin using fer- 
romagnetic electrical drain contacts that exhibit a spin- 
sensitive conductance [16]. Figure 7 demonstrates the 
utility of Eqs. (15)-(18) by computing all three compo- 
nents of electron spin polarization (n^, rix, and Uy) for 
the case of optically injected electrons (no || z) diffusing 
in the presence of both a real uniform magnetic field By 
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(left column), and in the presence of an effective chiral 
magnetic field due to strain (middle column). The 
rightmost column shows horizontal line-cuts through the 
center of the corresponding images of and rix. The 
diagrams at the bottom of Fig. 7 illustrate how elec- 
tron spins, generated at the center of the image, precess 
from initially out-of-plane (no || z) to an in-plane direc- 
tion as they diffuse. In the case of a uniform applied 
magnetic field By, all electrons precess in the same di- 
rection about By regardless of their momentum, so that 
Ux is finite with circular symmetry, and Uy is everywhere 
zero. In contrast, electrons diffusing in the presence of 
strain precess about the chiral effective magnetic field Be 
such that their spins are oriented radially away from the 
point of injection (at some characteristic radial distance). 
This radial spin distribution is reflected in the asymmet- 
ric images of Ux and Uy. In the presence of strain, the 
fact that the sign of electron spin polarization (for a spin 
component orthogonal to the spin generation term) is de- 
pendent on the direction of spin flow has been shown to 
be a valuable diagnostic tool (see Ref. [16]). 
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FIG. 7. 40 X 40 /xm images of the calculated z-, x-, 
and y-components of optically-injected spin polarization 
{nz,nx,ny). Spins are injected with no || z at the center 
of the images, and subject to diffusion only (no electric field) 
in the presence of either an applied magnetic field (left col- 
umn; By = 20 G), or strain (middle column; Cxy — 7 x 10~^). 
D = 10 cm^/s and Ts — 125 ns. The rightmost column shows 
horizontal line-cuts through the center of the corresponding 
images of riz and rix- The diagrams at the bottom illus- 
trate how spins precess from initially out-of-plane (no || z) 
to in-plane while diffusing in the uniform magnetic field By, 
or in the chiral effective magnetic field Be. 



For electrical spin injection from ferromagnetic con- 
tacts, the generation source is usually spatially extended 
and the orientation of injected spins follows the contact 
magnetization, which is typically in the x-y plane of the 
sample surface. However, polar magneto-optical Kerr mi- 
croscopy is sensitive to the z-component of spin po- 
larization. Therefore, a magnetic or strain field must be 
used to rotate the injected spin polarization out-of-plane 
in order to allow detection. In Fig. 8, we show calcu- 
lated images of resulting from a (theoretical) 40 jiia 
wide ferromagnetic stripe contact having in-plane mag- 
netization (M|| x). This spin polarization source is anal- 
ogous to the Fe/GaAs spin injection contacts used exper- 
imentally in Ref. [16]. The left column shows 100 x 200 
/im images of n^, and the right column shows horizon- 
tal line-cuts through the images. In Fig. 8(a), electrons 
injected with initial spin orientation no || x diffuse in 
the presence of strain {cxy = 1 x 10~^). Spins diffusing 
to the left and right precess into-plane and out-of-plane, 
respectively, giving the asymmetric image shown. The 
maxima and minima of nearly coincide with the edges 
of the source contact because the contact width greatly 
exceeds the characteristic spin diffusion length \/Dts. In 
Fig. 8(b), we consider the injected spins diffusing in the 
presence of a applied magnetic field {By = 1 G). The 
image of is symmetric, because electrons precess in 
the same direction (out-of-plane) regardless of which way 
they diffuse. Lastly, Fig. 8(c) includes an additional uni- 
form drift electric field {Ex = 4 V/cm), which shifts the 
spin polarization in the direction of electron drift (we 
consider the simplest case here, in which electric field is 
uniform across the image; i.e., the electric field is not 
generated by the injection contact, nor does the injec- 
tion contact distort the uniform electric field.) Note that 
finite spin polarization still exists on the "upstream" 
side of the source contact due to diffusion of injected 
electrons against the net electron current. 
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FIG. 8. 100 X 200 /xm images of the calculated z-component 
of electron spin polarization (uz) resulting from in-plane spin 
injection (no || M || x) by a 40-/im wide magnetized stripe 
contact (edges of the stripe are shown by dashed lines) . Corre- 
sponding horizontal line-cuts through the center of the images 
are shown to the right, (a) Spin diffusion alone for the case 
of a strained sample (cxy = 1 x lO"'^). Spins diffusing to the 
left or right precess in opposite directions, (b) Spin diffusion 
alone for the case of an applied magnetic field {By = 1 G). 

(c) Spin diffusion and drift for the case of applied magnetic 
and lateral electric fields (By = 1 G, Ex = 4 V/cm). For all 
images, /i — 3000 cm^/Vs, D = 10 cm^/s, Ts - 125 ns. 

Finally, we demonstrate the utility of the spin drift- 
diffusion equations in modeling the two-dimensional spin 
flows resulting from nontrivially-shaped injection sources 
(having either in-plane or out-of-plane magnetization), 
such as might be realized in actual experiments. Figure 
9 shows calculated images of in a semiconductor epi- 
layer resulting from a uniform square 80 x 80 /im injection 
contact. In Fig. 9 (a) and (c), the contact magnetization 
(and injected spin polarization) is oriented out-of-plane 
(M II no II z). For comparison, the contact magnetiza- 
tion (and injected spin polarization) in Figs. 9 (b) and 

(d) is in-plane and along the x-axis (M || no || x). Spins 
drift to the right in an electric field of = 8 V/cm in 
all of the images. 

In actual experiments, it is often desirable to mea- 
sure the spin polarization at a specific location on 
the sample as a function of in-plane magnetic field By 
[11,16]. These "Hanle curves" [40] provide valuable in- 
formation on the dynamics of the flowing spins, and in 
certain cases can be used to infer spin drift velocities and 
also, in conjunction with controlled strain, the direction 
of spin flow [16]. Figures 9(e) and (f) show calculated 
Hanle curves acquired at 4, 20, 40, and 80 fim from the 
right edge of the square injection contact. Solid (dashed) 
lines show the Hanle curves in the absence (presence) of 



strain. For the case of out-of-plane injection [Fig. 9 (e)], 
the Hanle curves are symmetric with respect to By. In 
contrast, for the case of in-plane injection [Fig. 9(f)], the 
Hanle curves are antisymmetric with respect to By , as ex- 
pected for orthogonal spin generation and detection axes. 
Unlike conventional Hanle measurements which are typi- 
cally sensitive to the spin polarization of all injected spins 
(e.g., the case of optical spin injection and detection via 
spatially- integrated luminescence polarization), in these 
studies the injection source and the point of detection 
may be spatially separated, so that only a subset of the 
injected electrons are detected by the probe laser. Fur- 
ther, these experiments allow access to a regime wherein 
the spin drift length exceeds the characteristic spin diffu- 
sion length. In this regime the Hanle curves can exhibit 
highly atypical shapes, even showing multiple oscillations 
as By is swept [16]. For the images of Fig. 9, where the 
injected spins subsequently flow to the right, the effect 
of strain is to shift the Hanle curves to negative values of 
By (dashed lines), implying that the effective magnetic 
field due to strain (Be) is oriented along -\-y. Conversely, 
a shifted Hanle curve can be used to infer the presence of 
strain in the sample, and also can be used (in conjunc- 
tion with controlled strain) to infer the actual direction 
in which the spin polarized electrons are flowing [16]. 
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FIG. 9. Calculated images of riz resulting from either 
out-of-plane (left column) or in-plane (right column) spin in- 
jection and subject to either magnetic or strain fields. The 
300 X 200 /xm images show the z-component of calculated 
electron spin polarization (uz) resulting from a square 80 x 
80 /im injection contact (edges of the source are indicated by 
dashed lines). In (a) and (c) the contact magnetization is 
out-of-plane (no || i), and in (b) and (d) the contact mag- 
netization is in-plane (no || x). Measuring Uz at a specific 
position as a function of in-plane magnetic field By yields the 
"Hanle curves" [40] shown in (e) and (f). These curves are 
acquired at distances of 4, 20, 40 and 80 //m from the right 
edge of the square source contact. Solid lines correspond to 
sweeping By alone; dashed lines include the presence of strain 
{exy = 0.03%), which shifts the Hanle curves to the left. In 
the calculations, fi = 3000 cm^/Vs, Ea: = S V/cm, D = 3 
cm^/s, and Ts = 125 ns. 



V. SUMMARY AND CONCLUSION 



We investigated the effects of electric fields, magnetic 
fields, and off-diagonal strain on the transport of spin- 
polarized electrons in zincblende semiconductors. Start- 
ing with a quantum kinetic approach, we first derived a 
semiclassical kinetic theory of electron spin dynamics and 
transport, and from this kinetic theory constructed spin- 
drift-diffusion equations when the total electron density 
is position independent and for spatially uniform electric, 
magnetic, and strain fields. This case of spatially uni- 
form fields and uniform electron density corresponds to 
the experimental situation realized in Refs. [15] and [16]. 
We compared the results of our spin-drift-diffusion model 
with Kerr microscopy images and found very good agree- 
ment. Our semiclassical kinetic theory is formulated on a 
microscopic foundation and can be extended to the cases 
of spatially varying (electric, magnetic and strain) fields 
and electron density. We contrasted the spin precession 
resulting from magnetic and strain fields. We found that 
because spin-strain coupling depends linearly on electron 
wave vector whereas spin-magnetic field coupling is inde- 
pendent of electron wave vector k, spatial coherence of 
electron spin precession is much better maintained with 
strain than with magnetic fields. This effect is most dra- 
matic for 1-D systems. In contrast to the case of a mag- 
netic field, we find that the spatial period of electron spin 
precession is independent of the applied electrical bias in 
samples with off-diagonal strain. The freedom to oper- 
ate at variable electric bias may benefit functional de- 
vices that are based on rotation of spin from one point in 
space to the other, as in spin transistor designs. Finally, 
we explore the utility of the spin drift- diffusion model by 
considering spin flows in a variety of potential cases, such 
as in-plane and out-of-plane spin injection from spatially 
extended sources in the presence of electric, magnetic, 
and/or strain fields. 
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VII. APPENDIX: SEMICLASSICAL KINETIC 
THEORY OF ELECTRON SPIN DYNAMICS 
WITH ELECTRIC, MAGNETIC AND STRAIN 
FIELDS 

We start from the Keldysh Green's function descrip- 
tion [42] for a system driven out of equilibrium [43] . The 
Dyson equation 



(Go^-E)0G = ^(1,2) 



(23) 



is written for the 4x4 matrix Green's function in the ten- 
sor product of spin and Keldysh space 



G 



K 



G« G 
G^ 



(24) 



where G^, G^ and are the retarded, advanced and 
Keldysh Green's function which are 2x2 matrices in spin 
space (for more details on this Larkin-Ovchinnikov nota- 
tion see Ref. [44]). The operator Gq ^ is defined as 







[ift^t, -e(lMl,2), 



(25) 



where the indices 1,2 represent the first and the second 
set of coordinate, time and spin variables and 5 {1^2) rep- 
resents a product of delta- functions in these variables; the 
operation (g) denotes matrix product in relevant spaces. 
Here £ is the self-energy operator while e expresses the 
effective mass Hamiltonian for the two spin states near 
the conduction band minimum. It is diagonal in Keldysh 
space and given by [40] 

e 0\ 

e ; ' 



^j2mFE, 



zK, + a{k X n) . 



The last term in fl^-^-^ with the coefficient a is the 
Rashba term (for an asymmetric heterostructure with 
symmetry-breaking axis n), while 



(26) 
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as 



= [Gp-i - S ? G]_ (T, R; cu, p) = i[G^^ - t, G]^ - 
[|0«^/(R; p) ■ <T © i + t{T, R; lo, p), G(T, R; lo, p)]_, (31) 



where 



[M, = -{d^Md^N - d^Md^N) 



is the Dresselhaus term arising from bulk inversion asym- Using the gradient expansion Eq. (29), Eq. (27) is written 
metry, ac = (4/3)(m/mc7;)?7/ yl — r//3, rricv determines 
an interband term quadratic in k and is related to the in- 
teraction of the conduction and valence bands with other 
bands. The components and are obtained by cyclic 
permutations of Eq. (26). From symmetry considera- 
tions alone [40,45] ^^^^ could also contain terms that 
couple spin to diagonal strain components, e.g., terms of 
the form = kz{exx — ^yy)^ ^tc. However, these terms 
only appear if spin-orbit mixing of states widely sepa- 
rated in energy is included [see Ref. [40] for more details] 
and therefore their numerical coefficients are small. 

The Dyson equation gives an exact description of the 
system. However, it cannot be reduced to an expression 
involving equal- time Green's functions which determine 
physical quantities, such as densities and particle cur- 
rents, via the single-particle density matrix. Deriving the 
quantum-kinetic equation involves considering the differ- 
ence of the Dyson equation and its conjugate. 



(32) 



is the generalized Poisson bracket and denotes a tensor 
product. 

Within the gradient approximation, we introduce 
gauge- invariant quantities as 



(33) 



[Go ^ - S f G]_ = 



(27) 



and approximations including a gradient expansion and 
a quasiclassical approximation. (Here [A^ B]- denotes 
a commutator of quantities A and B with respect to the 
operation (g).) We use a result due to Langreth [46] within 
the gradient expansion to obtain the kinetic equation for 
a gauge invariant distribution function. 

To introduce the gradient expansion, we transform to 
a mixed Wigner representation: 

R = ii±ia, T = r = n - r2, t = ti - t2 . 

After performing a Fourier transform with respect to rel- 
ative coordinates r and t 



A{n,F,TL,f)=A{Lj,p,-R,T) 
by a change of variables 

P = p - — (R, T), n = w - (e/h)4>{Ti, T) . (34) 
c 

The relations between the derivatives in old and new co- 
ordinates are 



edA 



dR-i 



G{X,p) 



dre-'P-''G{X^x/2,X-x/2) , (28) 



where X = (T, R), x = (t, r), p = (cj,p//i) and the dot 
product is given hj p-x = —ujt -h p • r//i, we represent the 
Fourier transform as a Taylor expansion expression in- 
volving products of those quantities in the Fourier trans- 
formed mixed representation: 

{M^N){X,p) = e'^^^<-^^-^^^^^M{X,p)N{X,p) . 

(29) 

We defined the exclusive partial derivatives dx = 
(— 9t,Vr), = {—duj^hdp) acting only on the quan- 
tity M; the scalar product of derivatives amounts to 
dx ' dp = —drduj + hdjidp. In the mixed representation 
the operator Gq^ is given by 

Go '(R; io,p)=hLo-^{p- -Af - e.^(R) 



and 9p = 9p, d^j = Oq^ where E is the electric field. 
From these relations one can prove the identity connect- 
ing the generalized Poisson brackets in old and new co- 
ordinates [46] 

[M,7V]^ = [M,7V]p + 
f {[B . (V^M X V^TV) + E . {d§M\/§N - \/§Md§N)] 
- [N ^ M]} . 

The generalized Poisson bracket in the new variables is 
given by 



^[Go\G]p = dTG + vdnG 
f E.J M[dpA"^9R,GU - [dR^nt",dp,G]+), 



(35) 



(30) 



where v = P/m and [A, 5]+ denotes an anticommutator. 
The Dyson equation in gauge invariant form becomes 

{dT + vdn + eE-{dp + ^dn)}G 
+ \Y.i{W^^p^f^ ,d^G\+ - [a.d^nff ,dpG] + ) 
+ £(v^// x'B)-&pG + {[n^ff •o-,G]_ =-i[S,G]_, (36) 

where v'^ =\—{h/2)aidp^1^^ ■ Equation (36) is a cen- 
tral result of our microscopic derivation. It describes spin 
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dynamics in the presence of electric, magnetic and strain 
fields provided the lengthscale of spatial variation is large 
compared to the characteristic lengthscale for wavefunc- 
tion variation. For the case of uniform strain, the strain 
terms in Eq. (36) agree with those found in Ref. [35]. 
In the absence of strain fields, the half-sum of the A 
and K components of Eq. (36) in the equal-time limit 
yields the kinetic equation obtained in Ref. [27] for the 
distribution function N = -(i/2)(G^ + - G^) [43]. 

The diagonal components G^ and G^ of the Green's 
function characterize the states, while the Keldysh com- 
ponent contains information on the occupation of 
these states [44]. We consider the Keldysh component 
of Eq. (36) in order to derive the kinetic equation. A 
distribution function /(T, R;0, P) is introduced by the 
ansatz 



(37) 



Using [Go ^-E^ ? G^]_ = and [Gq ^-E^ ? G^]_ = 0, 
the Keldysh component of Eq. (36) gives 



G^(8)B-B(8)G^ = 0, 



(38) 



where the 2x2 matrix B is given by 

B[f] = [Go ' f /]_ - ® / + / + 

^[Go-^?/]--C(S,/), (39) 

and in the mixed representation we have 

+ {dT+v-dR + eE-{dp + ^dn)}f, (40) 

where the Lorentz force term was neglected. In the geom- 
etry that we consider (Fig. 1), electrons are confined to 
move in a two-dimensional layer perpendicular to [001], 
and electric and magnetic fields are applied in the plane 
of electron motion. Hence the Lorentz force is along the 
z-axis and does not affect the motion in this case. 

We seek the solution of Eq. (38) for the distribution 
function as the solution of the equation B[f] = i.e. 



[Go-^?/]-=C(S,/) . 



(41) 



Taking the equal-time limit (which amounts to integrat- 
ing over Q so that terms with derivative vanish), 
Eq. (41) is then written as 



{dr -\- V • dn -\- eE • dp}h 

+ ^[n^ff . a, = - i / dnC{±, /), (42) 
where h(T, R, P) = / dnf{T, R; 17, P). 



For the experimental situation that we consider the 
applied electric, magnetic, and strain fields are spatially 
uniform and the net electron density is constrained by 
electrostatics to be spatially uniform. It is convenient 
to write h = po -\- Yl^=i Pi^i- For the steady state, spa- 
tially uniform case considered experimentally, Eq. (42) 
for the spin components of p becomes a Boltzmann type 
equation of the form 



(43) 



where Pc is a collision integral representing the right 
hand side of Eq. (42). To construct a spin-drift-diffusion 
equation from Eq. (43), we make the ansatz 



p,(R, P) = n,(R)/0(P) + ^,(R, P), 



(44) 



where ni(R) is a spin density that is to be determined, 
/^(P) is a known equilibrium momentum distribution 
function and ^^(R, P) is a small correction. The collision 
integral is described using a relaxation time approxima- 
tion 



(5i(R,P) 



(45) 



where r is a momentum scattering time that may depend 
on P but not on R. Substituting Eq. (44) into Eq. (43) 
and neglecting contributions for the correction term 5i 
on the left hand-side gives a first approximation for 6i 



= -r 



P dnA ,0 f ^ df\ 



(46) 



with cyclic permutation of indices for the other compo- 
nents of 6i. This procedure can be iterated to give a sys- 
tematic expansion for Si , that is Si = -\-Sj -\- ... with the 
momentum scattering time r being the small parameter 
in the expansion. We consider the case of rapid momen- 
tum scattering and keep only the lowest order term. 
We integrate Eq. (43) over all momenta to give 



9] 



- J d^P[n'^^ xp]i = J d^PiPc) = {Gi - Ri). (47) 

The second term on the left-hand side of Eq. (47) van- 
ishes identically and the possibility of spin generation 
(Gi) and relaxation {Ri) has been included empirically 
in the momentum integration of the scattering integral. 
[A more rigorous discussion of how momentum scatter- 
ing (Elliot- Yafet mechanism) and the Dresselhaus term 
in the Hamiltonian (Dyakonov-Perel mechanism) lead to 
spin relaxation within this formal structure is given in 
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Ref. [27].] Substituting Eq. (44) into Eq. (47) and using 
Eq. (46) for gives 

+ J d^PT^lff [(2^ • dnn,) f + n,(eE • dpf)] 
- /d^P^Oe// [(2£ . dnny) f + UyieE ■ ch>f )] 

= G,-^ (48) 

with cyclic permutations of indices for the other compo- 
nents of rii. Here the spin relaxation was taken to have 
the form Ri = n^/r^ where Tg is a spin relaxation time. 
Using the form given in Eq. (1) for ft^-^-^ gives Eqs. (4)- 
(6) with the mobility and diffusivity given by 

M = ey"d3pr£^((9e/°) (49) 

and 

where de denotes a derivative with respect to energy. The 
first term on the left-hand side (LHS) of Eq. (48) gives the 
usual diffusion and drift terms. The second term on the 
LHS, which does not contain the momentum scattering 
time T, is even in momentum for the magnetic field con- 
tributions to ft^-^-^ and odd in momentum for the strain 
contributions. Thus, the magnetic field contribution to 
this term survives the momentum integration but the 
strain contribution does not. By contrast the third and 
fourth terms on the LHS contribute for the strain terms 
but not for the magnetic field terms. The fifth term on 
the LHS is even in momentum for both magnetic field 
and strain terms, but odd for the cross term. For the 
strain terms there are quadratic-in-momentum contribu- 
tions from ft^-^-^ that combine with the momentum scat- 
tering time T to give a diffusivity whereas the magnetic 
field terms are proportional to the averaged momentum 
scattering time. Because the momentum scattering time 
is short, these magnetic field dependent terms are very 
small and can be neglected. For the experimental geom- 
etry ni^J^ is zero and ^%^^ /Py = fll^^/P^. Therefore, 
the contribution of this term in the equation for riz is 
twice that in the corresponding equations for Ux and Uy. 
This accounts for the factor of two multiplying C^D in 
Eq. (6) compared to Eqs, (4) and (5). For strictly 1-D 
motion in the x-direction ^y-^-^ vanishes and the factor 
of two multiplying C^D in Eq. (6) is replaced by unity. 
The sixth term on the LHS vanishes for the geometry 
that we consider. 
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